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Abstract 

The problem of sheath formation in front of a conductive planar plate in- 
serted into the plasma is formulated. Initially, the plate is assumed to be 
neutral. It is shown that the charging-up process of the plate is accompanied 
by the excitation of electron plasma waves. 

1 Introduction 

Investigations of sheath formation in front of a floating plate have hitherto been 
restricted to fluid studies on the ion time scale [1]. By contrast, the response of the 
plasma in the very early stages of sheath formation is not well known. In this paper, 
we present PIC simulations of the plasma dynamics over just a few electron plasma 
periods after the beginning of the process. These simulations have been performed 
by means of the BIT1 code [2], developed on the basis of the XPDP1 code from U. 
C. Berkeley [3]. 

A floating plate is placed in contact with a uniform, quasi-neutral plasma, which 
is assumed to be infinitely extended on one side. Due to the higher thermal velocity of 
the electrons, the plate starts charging up negatively, so that electrons are gradually 
repelled, ions are attracted, and a positive-space-charge sheath begins to form. An 
electron plasma wave is observed the properties of which strongly depend on the 
plasma characteristics (electron and ion temperatures, plasma density, etc.). 



Our PIC simulations are performed with different numerical set-ups and plasma 
characteristics. A full set of simulation diagnostics is used to measure the properties 
of the electron waves. 

2 Simulation set-up 

We consider a one-dimensional system. The planar conducting plate and the 
(artificial) right-hand boundary of the systems are placed at positions x = and 
x = L > 0, respectively. The length L is to be chosen large enough for this system 
to reasonably approximate a semi-infinite plasma (L » Xd, with A^ the electron 
Debye length). In order to have adequate resolution in space, the length of the grid 
cells has been selected as Ax ps Ad/2. 

2.1 Initial and boundary conditions 

At the initial time (t = 0) the electron and ion densities are equal (n e0 = n i0 = 
no), the distribution functions of both particle species are fully Maxwellian, and the 
electric potential is zero (V = 0) everywhere in the system, including the plate surface. 

Throughout the entire simulation, the following boundary conditions are applied 
to the particles: At the plate, all particles impinging are absorbed and no particles 
are injected into the plasma. At the right-hand boundary, on the other hand, all 
particles impinging are absorbed but new particles with half Maxwellian distribution 
functions are injected at a constant rate. The system is floating, i.e., the sum of 
particle plus displacement currents equals zero. According to these conditions we 
observe the following behavior. 

In the unperturbed plasma region (i.e., for x — > L) the electron velocity distri- 
bution function will not change appreciably (so that E(L,t) = 0), whereas at the 
plate it will acquire a cut-off form. This is because the negative-velocity electrons 
are absorbed by the plate and charge it negatively; during this process, the ions can 
be considered to be at rest. With increasing negative surface charge, the negative 
potential drop in the region close to the plate becomes higher and more and more 
electrons are reflected towards the plasma. After some time this perturbation propa- 
gates into the system. The shape of the distribution function essentially depends on 
the potential drop at the plate. 

Due to the loss of particles by absorption at the plate, the total number of particles 
in the system is dropping all the time. However, this aspect is not of great concern here 
because the total loss of particles during the entire simulation presented is negligible. 

2.2 Simulation parameters 

In the following tables we present the parameters used for our simulation. The 
(electron and ion) particle fluxes corresponding to the unperturbed plasma region 
are: 

T s = -^=-n 0s s = e,i (1) 
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These expressions are used to calculate the particle injection fluxes from the right- 
hand boundary. 



Plasma parameters 


Parameter 


Value 


Remarks 


J-e 


U.o ev 




T 
J-i 


U.i ev 




n i = n 0e = n 


0.0 X 1U m 


at i = U 


Vth e 


z.yo4y x iu m/s 




D-n, 

u thi 


3 095 x 10 3 m/s 




r e 


1.006 x lO 20 ™' 2 ^ 1 






1.05 x lO 18 ™- 2 ^- 1 






1.64 x 10 9 s _1 


electron plasma frequency 




3.83 x 10 7 s" 1 


ion plasma frequency 


rrii 


1.67 x 10~ 27 Kg 


proton mass 


Ad 


1.8 x 10~ 4 m 





Simulations parameters 


Parameter 


Value 


Remarks 


Ax 


8 x 10~ 5 m 


grid-cell length « 


L 


0.16 m 


system lenght 


S 


10~ 4 m 2 


plate aria 


At 


5.398 x nr n s 


time step 




1.079 x 10" 7 s 


total simulation time 



3 Simulation results 



Figure Q shows the potential profile 
close to the plate at t = 6.477 x 10~ 8 
s. The potential drop at the beginning 
of the sheath evolution is monotonic in 
space. After quick acquisition of nega- 
tive charge, the plate repels the electrons 
in the form of a pulse leaving behind a 
positive-space charge region. As a result, 
the potential close to the plate becomes 
slightly positive. In front of this region, 
the negative space charge produced by 
the primary-pulse electrons leads to a 
potential minimum ("virtual cathode"), 
which gradually reflects more and more 
slower electrons back into the plasma. 
These latter electrons spend a long time 
in the region of the virtual cathode and 



t = 6-4776x10"' s 




Figure 1: Potential profile in the perturba- 
tion region at t = 6.4776 x 10 -8 s 
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hence deepen its potential further. According to Figures. El and El this first potential 
perturbation (consisting of a potential hill and a potential well) propagates into the 
unperturbed plasma, with additional similar wave structures forming behind it. 



9 ° 



0.010 _, 
0.005 — 
0.000 - 

-0.005 .: 

-0.010 - 

0.010 

0.005 

0.000 
-0.005 
-0.010 



t =Atco =8.9e-2 

pe 



1200 i 



\ ' i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 r 



1100 i 



0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.010 
0.005 
0.000 
-0.005 
-0.010 

0.01 
0.00 
-0.01 
-0.02 



\ < i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 i 1 r 



10001 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



900! 



\ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



800 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



700 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 p 



600 1 



\ < 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



500 : 



\ < 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



400 : 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' 1 1 1 1 1 1 r 



300 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




200! 



I ' I 1 I ' I 1 I 1 I ' I 1 I ' I 1 I 1 I ' I 1 I 



100 1 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 

20 40 60 80 100 120 140 160 180 200 220 240 



Figure 2: Electron waves evolution 
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t = 6.4776x10 " s 
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To verify that these waves are phys- 
ical and not just due to numerical ef- 
fects, we have performed other simula- 
tions with different parameters. In par- 
ticular, we concentrated on the electron 
temperature. We know that the Debye 
length is proportional to the square root 
of the electron temperature. Hence, if we 
increase the temperature by a factor of 
four, the Debye length must increase by 
a factor of two. Since, in addition, there 
is a relation between the wavelength of 
the electron waves and the Debye length, 
the variation of the electron temperature 
should also have an effect on the wave- 
length. This is clearly illustrated in 

Figure where the wavelength is seen to increase with the square root of the electron 
temperature. 



Figure 3: Electric field profile (perturbation 
region); T e = 0.5eV 




Figure 4: Electron waves in plasma with different temperatures for electrons T e = 2eV 
and T e = 0.5eV 



4 Summary and conclusions 

This work represents the beginning of a self-consistent kinetic study of sheath 
formation, taking into account both electron and ion dynamics. Here, during the 
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short simulation time considered, the ions are practically immobile, and only the 
electrons take part in the process. In the next step, the effect of ion dynamics on 
sheath formation will be considered as well. 
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